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Abstract 



There are a number of mathematical theorems in the hterature 
on the dynamics of cosmological models with accelerated expansion 
driven by a positive cosmological constant A or a nonlinear scalar field 
ly-^ ■ with potential V (quintessence) which do not assume homogeneity and 

. isotropy from the beginning. The aim of this paper is to generalize 

CJ ! these results to the case of fc-essence models which are defined by a 

Lagrangian having a nonlinear dependence on the kinetic energy. In 
^H^' particular, Lagrangians are included where late time acceleration is 

driven by the kinetic energy, an effect which is qualitatively different 
from anything seen in quintessence models. A general criterion for 
^ ■ isotropization is derived and used to strengthen known results in the 

^ . case of quintessence. 



1 Introduction 

In recent years many models have been introduced as possible explanations 
of the observed acceleration of the expansion of the universe. The simplest 
and best-known of these involve a positive cosmological constant A or a non- 
linear scalar field with potential V. There is a huge literature on this subject 
- see for instance for a review. In most papers the analysis is restricted 
to homogeneous and isotropic cosmological models since it is expected that 
these will give a good approximate description of more general spacetimes. 
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For some of the simplest cases with late time accelerated expansion this ex- 
pectation has been confirmed by rigorous proofs. There are a number of 
papers where mathematical results on the dynamical properties of models of 
this kind are obtained in the homogeneous and inhomogeneous cases. Ho- 
mogeneous models with A > were studied in |2ZI and [T7j and those results 
were extended to scalar fields with potentials belonging to various classes in 
[T^ , PH] 5 |2I] and |23] ■ In the inhomogeneous case theorems were proved for 
vacuum solutions with A > in [9j, |j221 and PP, the case of an exponential 
potential in |T3| and for some inhomogeneous cases with matter in jJE] and 
Note also the results of on curvature coupled scalar fields. The 
purpose of this paper is to extend some of these results to a larger class of 
models, the fc-essence models. 

The name fc-essence denotes a nonlinear scalar field (f) with Lagrangian 
where X = — |Va0V"0, which is used to describe late-time cosmic 
acceleration. It was introduced by Armendariz-Picon, Mukhanov and Stein- 
hardt 0] following the earlier application of similar models to the very early 
universe under the name fc-inflation j2] . Fields of this kind arise in effective 
field theories coming from string theory when higher order corrections in the 
string tension or number of loops are included. Special cases of the /c-essence 
Lagrangian are the ordinary minimally coupled nonlinear scalar field (asso- 
ciated with the name quintessence) where L{(j),X) = —V{(f)) + X and the 
tachyon, where L = —V{(t>)^/1 — 2X (see |Hj, [II])- Here the potential ^ is a 
smooth function, assumed in the following to be non-negative. The equation 
of motion is 

In the spatially homogeneous case this reduces to 

Note that in the spatially homogeneous case X > 0. To ensure the solvability 
of equation for it should be assumed that dL/dX + 2Xd'^L/dX'^ ^ 0. 
For a given choice of L this is a restriction on the allowed values of and (f). 
In general L is not defined on the whole (0, X)-plane but only on a subset 
of it. As an example, for the tachyon it is necessary that X < 1/2 in order 
that L be defined. Once this is ensured dL/dX + IXd'^LjdX'^ > under 
the sole restriction that the potential be positive. 
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The energy-momentum tensor is given by 
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T"^ = —V"(t)V^(t) + Lg''^ (3) 
oX 

In the spatially homogeneous case it has two distinct non-vanishing compo- 
nents, p and p. In fact the pressure p is just equal to the Lagrangian. The 
energy density p is given by 2XdL/dX—L. It follows that p+p = 2XdL/dX . 
The inequality p + p > is a not very restrictive energy condition. In partic- 
ular it follows from the dominant energy condition. If this energy condition 
is satisfied then it follows that dL/dX > in the homogeneous case. In the 
following it will always be assumed that 

dL ^a'^L dL 

As is discussed in Section 2 these conditions follow from the hyperbolicity of 
the equation of motion of and the energy condition mentioned above. The 
sign of p + 3p = 2{XdL/dX + L) determines whether there is accelerated 
expansion in the homogeneous case. It follows from the equation of motion 
that 

(«) 

The paper is organized as follows. Section |21 deals with the question, for 
which choices of the Lagrangian L the equation of motion of fc-essence is 
hyperbolic and under what conditions on L various energy conditions are 
satisfied. In Section IHl late time acceleration and isotropization are shown for 
homogeneous spacetimes of Bianchi type I- VIII with fc-essence defined by a 
Lagrangian chosen from a wide class. The conditions defining this class are 
modelled on those previously used in analysing nonlinear scalar fields with a 
potential which have a positive lower bound. In particular tachyon models 
whose potential has a positive lower bound are treated. In these results the 
spacetime is allowed to contain normal matter in addition to the fc-essence. 
The normal matter is only supposed to satisfy the dominant and strong 
energy conditions. Section 0] proceeds to consider homogeneous spacetimes 
where late time accelerated expansion is driven by the kinetic energy. In these 
spacetimes the kinetic energy does not go to zero at late times, in contrast to 
those studied in Section El In Section a general criterion for isotropization 
is obtained and is applied in the special case of quintessence to generalize 
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the results of j2Sl- Appendix 1 contains some technical background required 
for Section 121 while Appendix 2 contains continuation criteria which are used 
to obtain global existence theorems in Section El and |31 

2 The general inhomogeneous case 

The equations of a fc-essence field coupled to the Einstein equations form 
a well-posed system. This means that if suitable initial data are prescribed 
there exists a unique local solution of the Einstein-matter equations inducing 
those initial data. This can be proved following the general scheme described 
in ^Uj, section 5.4. Under the assumptions (@)) the equation of motion of (f) 
is a nonlinear wave equation. This follows from the fact that the tensor 
multiplying the second derivatives of in (Q) is a Lorentz metric. The proof 
of this is given in Appendix 1. By reduction to first order this wave equation 
can be written as a symmetric hyperbolic system. The coupling of this system 
to the Einstein equations is such that the whole Einstein-matter system can 
be written in symmetric hyperbolic form. A local existence and uniqueness 
theorem follows. This argument also applies to the case where in addition 
to the coupling to the fc-essence field the Einstein equations are coupled to 
other matter having a well-posed initial value problem, e.g. a perfect fluid 
with reasonable equation of state. 

One cautionary remark is necessary. For the coupled system initial data 
may only be prescribed on a hypersurface which is spacelike with respect 
to both the spacetime metric and the metric in the principal part of the 
wave equation for 0. If d^L/dX"^ > then it is enough to assume that the 
hypersurface is spacelike with respect to g^/s- Otherwise it is not. In the 
latter case superluminal propagation of signals is possible. This is proved in 
Appendix 1. 

Energy conditions can be investigated on a case by case basis as is done 
in Appendix 1 for the signature of the metric in the equations of motion. 
The results dbTQ clS follows. The dominant energy condition is satisfied iff 
dL/dX > and XdL/dX — L > 0. For X > the weak energy condition 
is equivalent to dL/dX > and 2XdL/dX - L > 0. For X < the weak 
energy condition is equivalent to dL/dX > and L < 0. The strong energy 
condition is equivalent to dL/dX > and XdL/dX + L > for X > and 
dL/dX > and -XdL/dX + L > for X < 0. 

A potential difficulty with /c-essence models is that the matter may form 
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singularities which have nothing to do with the famihar spacetime singulari- 
ties of general relativity. (Cf. the discussion in jS].) It would be interesting to 
know for which choices of L this can be avoided. A simple place to start is to 
look at a /c-essence field in Minkowski space and consider solutions evolving 
from data which are close to zero. A criterion for the corresponding smooth 
solution to exist globally in time is the null condition ^Hl- This is an alge- 
braic condition on the nonlinearity. It turns out that it is satisfied for the 
fc-essence field if and only if d'^L/dcf)^ = at {(p,X) = (0,0) for k = 1,2,3. 
In that case there are no singularities for small data. In particular this is 
true if L only depends on X. 

3 Solutions with a positive lower bound for p 

Consider a spatially homogeneous solution of Bianchi type I- VIII of the Ein- 
stein equations with matter. Suppose that the energy-momentum tensor is a 
sum of two parts. One of these T*^ represents ordinary matter and satisfies 
the dominant and strong energy conditions. The other represents dark en- 
ergy and matter quantities derived from it will be given the subscript DE. 
It is assumed in the following that the dark energy satisfies the dominant 
energy condition. Some basic equations are as follows: 

-H' - fipDE + trSoE) - l^at^"' - f{p'' + trS'^) (6) 
fpDE + l{a^,^-'-R) + ^-^p'' (7) 
-^TcipDE + ^tTSDE) - ^cj^ta''' + ]:R - Mp'' + ltr5*0 (8) 

6 2b 6 

These are the direct generalizations of equations (4), (6) and (7) of fH] . 
Suppose that H is zero at some time ti. Then by (|7j) pde is zero at that 
time. By (jSJ it vanishes at all later times. Since the /c-essence field satisfies 
the dominant energy condition it follows that it makes no contribution to the 
energy-momentum tensor. Hence the geometry and normal matter satisfy the 
equations in the absence of fc-essence. By uniqueness for the full system of 
equations this must be true at all times and not just after time ti. It may 
be said that the /c-essence field is passive and these solutions are physically 
uninteresting. The results of [21] show that the normal matter vanishes and 



dH 
"dt 
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dH 
"dt 
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that spacetime is flat. For this reason it will be assumed in the following 
that H is everywhere positive. 

In the next theorem no specific choice of a dark energy model is made. 
It will be specialized to the case of /c-essence later. 

Theorem 1 Consider a solution of the Einstein equations of Bianchi type 
1-Vlll coupled to a matter model representing dark energy and other mat- 
ter satisfying the strong and dominant energy conditions. Let the energy- 
momentum tensor T^^ of the dark energy have the algebraic form of that 
of a perfect fluid and satisfy the dominant energy condition. Suppose that 
the energy density of dark matter p^E has a positive lower bound po > 
and that the solution is initially expanding {H > 0). Then H > Hq for 
some constant Hq > and if the solution exists globally in the future the 
quantities , aahO^""^ and R decay exponentially. 

Proof The first conclusion follows immediately from the Hamiltonian con- 
straint. To prove the remainder of the theorem let Z = dH"^ — 2A7cp£)E- As in 
the special case of quintessence considered in j^I] the quantity Z satisfies the 
inequality dZ/dt < —2HZ. The only properties of the dark energy model 
which are used in deriving this are those included in the hypotheses of the 
theorem. From the inequality it follows that if the solution exists globally 
in the future the lower bound for H implies that Z decays exponentially as 
t — s> oo. The Hamiltonian constraint shows that 

Z = -^R + 247rp^ (9) 

and this concludes the proof. 

The interpretation of this theorem is that global solutions isotropize and 
that the spatial curvature and the energy density of normal matter have 
very little effect on the dynamics at late times. In the case of /c-essence the 
energy-momentum tensor of a homogeneous model always has perfect fluid 
form. To obtain the other conditions it suffices to supplement by the 
inequality 2XdL/dX — L > C for a constant C > 0. The theorem also 
applies to fluid models for dark energy which satisfy the dominant energy 
condition and have a positive lower bound for the energy density, such as the 
Chaplygin gas il4j. Note that the argument of the theorem also applies to 
the situation where there is no global positive lower bound for p^E but for 
a given solution which exists globally in the future pnE tends to a positive 
limit as t ^ cxD. 

It follows from equation (0) and the lower bound for H that the quantity 
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XdL/dX is integrable on the positive time axis for any global solution. Now 



d/dt X 



dL_ 

dX 



) 




dX^ 



) 



+ X 



d(j)dX 



(10) 



When the expression for (f) is substituted into this the contribution containing 



H is non-positive. Discarding it gives an upper bound for the time derivative 
of XdL/dX in terms of various other quantities. 

Now a list of relevant assumptions on the Lagrangian will be collected. 
Given a constant Ci > there exists a constant C2 > such that: 

1. if p < Ci then X < C2 

2. if p < Ci then \dL/dX + Xd'^L/dX'^\ < C2 

3. if X < Ci then \dL/dX + 2Xd^L/dX^\ > 

4. if p < Ci then \d^L/d(l)dX\ < C2 

5. if p < Ci then \dL/d(j)\ < C2 

6. if X < Ci then dL/dX > 62^ 

7. there are positive constants C3 and C4 such that X < C3 implies 
L((/.,X) < -Ci. 

Theorem 2 Consider a solution of the Einstein equations of Bianchi type 
I- VIII coupled to a /c-essence model and other matter satisfying the strong 
and dominant energy conditions. Let the energy-momentum tensor T^^ of 
the /c-essence satisfy the conditions (@]) and the dominant energy condition 
and let there be a positive lower bound for the energy density. Suppose 
that the solution is initially expanding {H > 0). Then if the solution exists 
globally in the future and conditions 1.- 7. above are satisfied the expansion 
is accelerated at late times. 

Proof Equation shows that is bounded above. Then assumption 1. 
above implies that X is bounded. Assumptions 2.-5. imply that XdL/dX — >■ 
as t ^ 00. By assumption 6. we can conclude that X ^ as t ^ 00. 
Finally assumption 7. implies that the expansion is accelerated at late times. 

Since the assumptions 1.-7. look quite abstract, let us see what they 
mean for quintessence and for the tachyon. Consider first quintessence with 
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a non-negative potential. Conditions 1.- 4. obviously hold. Condition 5. 
holds provided V is bounded on any interval where V is. This assumption 
played an important role in f5T. Condition 6. is obvious. Condition 7. holds 
if V has a positive lower bound. This special case of Theorem 2 was proved 
in mi 

Next consider the tachyon. Since X is a priori bounded by one in this 
case condition 1. is obvious. In the case of the tachyon if p is bounded then 
so is V. If in addition V has a positive lower bound then X is bounded away 
from one. Then condition 2. follows. If V has a positive lower bound then 
condition 3. holds. For conditions 4. and 5. it should be assumed that V 
is bounded whenever V is. With a positive lower bound for V condition 6. 
holds. Condition 7. is automatic for the tachyon. Thus Theorem 2. implies 
the following result 

Theorem 3 Consider a solution of the Einstein equations of Bianchi type 
I- VIII coupled to a tachyon field and other matter satisfying the strong and 
dominant energy conditions. Suppose that the tachyon potential V has a 
positive lower bound, that V /V is bounded and that the solution is initially 
expanding [H > 0). Then if the solution exists globally in the future the 
expansion is accelerated at late times. 

For a perfect fluid or collisionless matter the global existence assumption 
in Theorem 3 holds automatically. For 

d/dti\og{l - 0')) = (1 - 0')(3if02 + V'/V^) (11) 

and this prevents |0| approaching one in finite time. This in turn ensures 
that (0, 0) remains in a compact subset of the domain of definition of L. It 
follows that Lemma 1 and Lemma 2 of Appendix 2 can be applied. 

If L depends only on X and is globally defined and if it satisfies condition 
1. then conditions 2.- 6. are automatic. Moreover condition 7. can be 
replaced by the condition L{0) < 0. In ^2j the case where L only depends on 
X plays an important role but L does not satisfy the positivity conditions 
Q everywhere. An alternative interpretation is to keep those conditions 
and say that L is only defined locally. An interesting generalization is given 
by L{(j),X) = Lo(0) + Li(0)X + L2(0)X2. This is considered in j2] under 
the condition Lq = (no potential). Theorem 2 applies if —Lq and Li are 
bounded below by a positive constant. 
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4 Solutions with a positive lower bound for 
the kinetic energy 

The original aim of A;-essence models was to produce a scenario where the 
accelerated expansion of the universe is driven by the kinetic energy of a 
scalar field. This is very different from the usual case of quintessence where 
the potential energy is the driving force. In the models studied in the previous 
section it was the latter mechanism which was at work. The present section 
is concerned with models closer to the original motivation for considering 
fc-essence. They have the property that they have late time accelerated 
expansion where the kinetic energy does not tend to zero. 

What can happen is illustrated by the following specific example. Let 
L = X"^ — 2X on the interval J = (1, oo). The conditions (0} are satisfied on 
J, as is the dominant energy condition. The equation of motion of implies 

X = -Q{'iX -1)-^{X -l)XH (12) 

In particular X is decreasing. Also H > (87r/3)(3X^ — 2X). It follows that 
if a solution of the Einstein-matter equations exists globally in the future 
then X ^ 1 as t ^ oo. The continuation criteria of Appendix 2 show that 
when the normal matter is described by a perfect fluid or coUisionless matter 
solutions exist globally. The mean curvature remains bounded away from 
zero at late times because pde does. The inequality dZ/dt < 2HZ then 
implies as in the last section that at late times the influence of the spatial 
curvature and the energy density of ordinary matter becomes negligible and 
the geometry isotropizes. Since XL' + L < for X G (1,4/3) it follows that 
all solutions are accelerated at late times. 

This example can easily be generalized. Consider a Lagrangian L = L{X) 
defined on an interval {Xq, Xi) with Xq > and Xi finite or infinite. Suppose 
further that (@)) holds together with the dominant energy condition. Finally, 
suppose that 2XL' — L is never zero, that limx^Xo(^ ~ ^o)^^L' exists and 
that limx^Xo L exists and is negative. The equation of motion is 

X = -6(L' + XL"Y^L'XH (13) 

Thus X is decreasing. From the Hamiltonian constraint it follows that > 
(87r/3)(2XL' — L). Hence under the given conditions if a solution exists 
globally in the future then X — > Xq as t — >■ oo. Provided L' + XL" is 
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bounded away from zero in a neighbourhood of Xq global existence theorems 
follow from the continuation criteria of Appendix 2. In this case it can be 
concluded as before that the geometry isotropizes and that there is late time 
acceleration. 

In the general case where the Lagrangian also depends on the equation 
of motion implies that 

(14) 

As a global simplifying assumption, suppose that the domain of definition of 
L is of the form Jxl where J is an open interval and / is the interval (Xq, oo) 
with Xo > 0. Assume as usual (H)) and the dominant energy condition. 
In the following a theorem is proved which concerns a situation where the 
dependence of L on does not disturb the behaviour seen for a Lagrangian 
only depending on X too much. Looking at (fT^ suggests that a smallness 
assumption should be made on dL/d(j) and d'^L/d(f)dX. It will be assumed 
that there is a positive constant a such that 

|0(9L/(90| + 2\(l)Xd'^L/d(l)dX\ < a{dL/dX) (15) 

This is motivated by the consideration of Lagrangians of the form L(0, X) = 
L{X)(j)~'^ as introduced in The time derivative of has constant sign 
and only solutions for which > will be considered. Since a solution which 
exists globally in the future satisfies > y/2XQ it follows that grows at 
least linearly as t — > oo so that (|TH|) implies decay of {dL/ d(j)) / {dL/ dX) . The 
following conditions are required in the following results. Given a constant 
Ci > there exists a constant C2 > such that 

1. if X > Xo + Ci then dL/dX > 

2. (X — Xo)^^(9L/(9X converges to a limiting function as X ^ Xq 

3. L converges to a negative limiting function Lq as X — > Xq 

4. if X < Xo + Ci then {2XdL/dX - L) > C^^ 

5. {dL/dX + Xd'^L/dX^) > C^^ for all X 

6. dhjdX < C2XI/2 for all X 
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Theorem 4 Consider a solution of the Einstein equations of Bianchi type 
I- VIII coupled to a fc-essence model and other matter satisfying the strong 
and dominant energy conditions. Let the energy-momentum tensor T^^ of 
the fc-essence satisfy the conditions (jH) and the dominant energy condition. 
Suppose that the inequality (fT3jl and conditions 1.-4. above are satisfied. If 
the solution is initially expanding and exists globally in the future then the 
expansion is accelerated at late times and the geometry isotropizes. 
Proof Under the assumptions of the theorem the following estimate is ob- 
tained: 

QX^H + 2^2X^1' 4^ - V2X'/'^ > X'/'^{6X'/'H - V2r'a) 
dX d(l)dX d(j) - dX^ v r ; 

(16) 

The function ^ tends to zero as t — > oo. Hence the second term in the 
bracket on the right hand side tends to zero. To get information about 
the first term, note that due to (jH) the energy density pde is an increasing 



function of X and that H > ySnpDE/S. Together with assumption 4. this 

implies a global positive lower bound for p£,E and for H. Hence X < for 
t sufficiently large. It follows using assumption 1. that X Xq as t oo. 
The remaining conclusions of the theorem follow using assumptions 2. and 
3. 

It is of interest to ask when global existence holds so that the theorem 
can be applied. First a lower bound for X will be obtained. At a given time 
either X < 2Xq ot X > 2Xq. If the former condition holds let t2 be the 
last time before ti at which X > 2Xo if such a time exists and let ^2 be the 
initial time otherwise. On this interval X and are bounded. It follows 
by condition 5. that 

|X| < CdL/dX (17) 

for a constant C. Since dL/dX is 0{X — Xq) as X — > Xq it can be concluded 
that X remains bounded away from zero on the interval [t2, ti] with a bound 
depending only on the initial data. It follows that for any solution X remains 
bounded away from Xq on any finite time interval. Using assumption 6 above 
and Gronwall's inequality it can be concluded that X is also bounded on 
any finite time interval. It follows that the solution cannot approach the 
boundary of the domain of definition of L in finite time. The results of 
Appendix 2 can be applied to obtain global existence statements. 
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The conditions 1.-6. do not apply to the case L{(f),X) = (f)~'^L{X) occur- 
ring in Next some properties of solutions for Lagrangians of this form 
will be examined. The equation of motion takes the form: 

(L' + 2XL")X + L'{QXH - 4720-^X^/2) ^ 2v^X^/V"^^ = (18) 

In this case it is possible for X to become zero, the condition for this being 
that 2XL' - L = (3/72)0X^/2^1'. In a FLRW model without ordinary 
matter this reduces to L = 2X{—6tt{L'Y + L'). Any solution of this equation 
leads to a solution of the Einstein equations coupled to fc-essence alone in 
which X is constant and the scale factor has power-law behaviour. Solutions 
of this type were considered in section 5 of and correspond to the k- 
attractors of |3]. If XL' + L < the expansion is accelerated. Eliminating L 
gives the inequality L' > I/Att. The dominant energy condition is equivalent 
to XL' — L > which gives the inequality L' > l/127r. 

It is clear that there are many /c-essence Lagrangians together with solu- 
tions of the type just discussed where X is constant. Note however that it is 
not guaranteed that they will exist for any given Lagrangian. In the case of 
the tachyon solutions of this kind are known explicitly [7]. 

By introducing a new time coordinate we can consider (|T8|) as an ordinary 
differential equation for X alone which does not involve 0. This makes it very 
easy to do a local stability analysis of the stationary solution in the class of 
spatially flat Friedmann models without normal matter. It turns out that 
the condition for stability is L"{XL' — L) > 37r(I/')2. Assuming the dominant 
energy condition this can only be satisfied if L" > which means that there is 
no superluminal propagation in inhomogeneous models for this Lagrangian. 
It is more difficult to prove something about the stability of these solutions 
within a class of anisotropic solutions. One reason this is difficult is because 
the accelerated expansion is only of power-law type and not faster. In [23] 
isotropization was proved for models whose expansion increases in time faster 
than any power of t corresponding to a potential which decays slower than 
any exponential. For general potentials which decay at a rate comparable 
to an exponential no result of this kind was obtained. The only case for 
which isotropization has been proved is that of an exact exponential ^H] • In 
the next section some partial results on isotropization are obtained which, in 
particular, considerably extend what is known in the case of quintessence. 
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5 A criterion for isotropization 



Consider a solution of Bianchi type I- VIII of the Einstein equations coupled to 
dark energy (satisfying the dominant energy condition) and ordinary matter 
(satisfying the dominant and strong energy conditions). A computation using 
the evolution equation for H and the Hamiltonian constraint shows that 



dH 
It 



> -3H' 



Air 

1 

3 



Pde 



(19) 



Suppose that limini t^oo{PDE — I^^Sde) / H"^ > Then it follows that for 
t sufficiently large dH/dt > —jSH^ for a constant j3 < 1. Hence H{t) > 
f3^^t^^ + 0(t~^). It follows that for t large H{t) > '~ft~^ for any 7 less than 
The constant 7 may be chosen larger than one. Putting this inequality 
into the relation dZ/dt < -2HZ shows that Z = 0(^-^7) and that Z/H"^ -> 
as t — i> 00. Thus a criterion for isotropization has been obtained. 

What does the criterion look like in the case of quintessence? There 



Pde - ^tiSDE ] /H"^ 



(20) 



so that we get the inequality liminf (87ry^/3if^) > 2/3. It will now be shown 
that if a 



hmsup^^oo(~^'/^) < 4y7r/3 then the criterion is satisfied. The 
starting point of the proof is provided by equation (11) and the inequality (13) 
of 



There is a number 0i such that —V'/V < a/2 + y47r/3 for all > 0i. 

Now distinguish between the cases <p/H > 1/712^ and ^/H < l/VU^. In 
the first case equation (11) of j2S] gives an upper bound for d/d(f){3H^/87iV) 
while in the second case equation (11) and the inequality (13) of that reference 
imply the bound (14) for d/d(j){3H'^/87!'V). Combining these two bounds 
shows that the following inequality holds in general: 
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on some interval for a positive constant Ci 
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(21) 
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(23) 
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The expression on the right hand side of this inequahty is strictly greater 
than a / \/487r. Increase 0i if necessary so that 



-2 



a 



+ 



a 



V48 



7r 



1 + C^i 1 



n -1 



V48 



7r 



(24) 



for (j) > (pi. This is possible due to the restriction on a which has been 
assumed. It follows that the criterion for isotropization holds. 

The result just proved extends the proof of isotropization and decay of 
spatial curvature and energy density of ordinary matter significantly in com- 
parison with what is proved in [23j. It applies for example to the potentials 
V{(f)) = Vo(log 0)^0" exp(— A0), the special case of equation (18) in j2SI with 
m = 1 provided A < AJn/S and to a number of other potentials listed in 



6 Conclusions and outlook 

In this paper results are obtained on the dynamics of solutions of the Ein- 
stein equations with /c-essence and normal matter. They include anisotropic 
spacetimes of Bianchi type I- VIII and not only FLRW models. In this set- 
ting late time isotropization is proved. They allow the normal matter to be 
anything which satisfies the strong and dominant energy conditions. Thus 
the results are not confined to fluids or multifiuids and apply, for instance, 
to coUisionless matter. All the arguments are rigorous and not dependent on 
heuristics. 

Before coming to the global results for homogeneous solutions, basic in- 
formation is presented concerning the initial value problem for the Einstein- 
matter equations with /c-essence in general, i.e. without symmetry assump- 
tions. A link is made with a concept in the theory of nonlinear hyperbolic 
equations, the null condition, and it would be desirable to establish further 
connections in the future. The paper ^ could be a good starting point. 

It is shown that in general if satisfies a suitable lower bound homo- 
geneous solutions of Bianchi type I- VIII which exist globally in the future 
become homogeneous with dimensionless measures of the spatial curvature 
and density of normal matter tending to zero exponentially. This works if a 
global lower bound is available for the given matter model or if it is avail- 
able for a particular solution. It is shown further in Section 3 that if certain 
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bounds are assumed on the Lagrangian it is possible to prove late-time accel- 
erated expansion. In this class of models X — as t oo. They very much 
resemble quintessence models with a lower bound for the potential. It would 
be interesting to prove a similar result for a class generalizing quintessence 
models with a potentials tending to zero at infinity as analysed in |23j . 

In Section 4 results are obtained on late time behaviour of solutions whose 
kinetic energy does not tend to zero as t —>■ oo. They achieve the aim of 
proving late time accelerated expansion and isotropization for models where 
these phenomena occur. The questions of behaviour at intermediate times 
and late time behaviour without acceleration are not addressed since they 
are not amenable to the techniques used here. If the expansion is not acceler- 
ated at late times then the behaviour is strongly dependent on Bianchi type 
and a unified treatment is difficult or impossible. As far as the behaviour at 
intermediate times is concerned it is difficult to even formulate interesting 
and precisely defined statements susceptible to rigorous proof. This is con- 
nected to the fact that the interesting statements often involve essentially 
quantitative elements while rigorous results tend to be qualitative in nature. 
These issues deserve a separate discussion. 

The first case considered in Section 4 is that where L is independent of (p. 
A more general theorem is then proved which shows that certain kinds of de- 
pendence on do not change the qualitative behaviour in comparison to the 
case where L only depends on X. Late time acceleration and isotropization 
are proved. Next the case L{(j),X) = (f)~'^L{X) is considered. This allows 
solutions with X = 0, the fc-attractors. The stability of these solutions 
within the class of isotropic spatially fiat models without ordinary matter 
is determined. Isotropization is not proved. The fact that the expansion of 
the candidate attractor is only of power-law type makes this difficult, as it 
does for quintessence models with this type of asymptotics. To make further 
progress a criterion for isotropization in dark energy models is developed in 
Section 5. It is shown that this criterion leads to a significant improvement 
in the results which can be obtained in the case of quintessence. Similar 
conclusions have not yet been obtained for more general fc-essence models. 

This paper contains new results on the mathematical properties of solu- 
tions of the Einstein equations coupled to dark energy and other matter and 
concentrates on establishing a basis for the mathematical study of /c-essence 
models. The fc-essence Lagrangians give rise to a considerable variety of be- 
haviour and the task of obtaining an overview of the possibilities is a challenge 
for the future. 
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Appendix 1. Analysis of characteristics 

Results similar to those of this appendix were obtained in The treatment 
here is intended to make the proofs more transparent by avoiding the use of 
series expansions. 

Consider the expression §^g"^ — f^V"0V^0 which occurs in the prin- 
cipal part of the fc-essence equations of motion. Evidently this is degenerate 
if dL/dX = 0. We now assume that dL/dX ^ and define 



h^f^ = (dL/dXy^—g"^ - T^^'"(p^^(p) = g""^ - KV^cpV^cp (25) 



where K = ^^fyff^. If 1 - iCVa^V^ is non-zero let h^is = gap + i^ksjUv 
Then by a direct computation is the inverse of h""^ . The condition 
1 - KVa(pV''(f) = is equivalent to gV^^V^^ = If = then this 
never happens. Otherwise the assumption that it does not happen puts a 
restriction on Vcr0V°"0. 

Next the signature of the metric h"^ will be investigated. Suppose first 
that the gradient of is timelike. Then V^(f) = AT" for a unit timelike vector 
T°. Then h'^^TaTp = — 1 — KA^ while h'^^SaSp = 1 for any unit spacelike 
vector S*" orthogonal to T'^. Moreover h'^^SaTp = 0. Thus the signature of 
h is Lorentzian for KA^ > —1. For KA^ < —1 it is positive definite. The 
equation becomes elliptic. If the gradient of (p is spacelike then it follows 
that = AS" for a unit spacelike vector S". A calculation similar to 
the one just done shows that the signature is Lorentzian for KA"^ < 1 and 
(— , — , +, +) for KA^ > 1. If V"0 is null the signature is always Lorentzian. 
Hence in all cases an equivalent condition to the metric being Lorentzian 
is ii'Va^V"^ < 1. This can be re-expressed as 2X^^ ~^ §k ^ ^ when 
dL/dX > 0. 

Let v" be a vector which is null with respect to the metric hajs- Then 
g^pv"v^ = -i^(V„0w")V(l - i^V„0V"0)^. It follows that if K is positive 
the null cone of haf3 is inside (or on) that of gap while if K is negative it is 
outside (or on) that of gap- 

Appendix 2. Continuation criteria 

In this appendix basic local existence theorems and continuation criteria 
will be proved for homogeneous spacetimes containing fc-essence and normal 
matter. The choices of normal matter discussed are a large class of perfect 
fiuids with positive pressure, a mixture of several non-interacting perfect 
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fluids and collisionless matter described by the Vlasov equation. Let G be 
the open subset of where the Lagrangian L of the fc-essence field is defined. 
Assume that the equation of state p = f{p) of any perfect fluid occurring is 
such that / : [0, oo) [0, oo) is continuous with /(O) = and that for p > 
the function / is with < /'(p) < 1 

Lemma 1 Consider a solution of Bianchi type I- VIII of the Einstein equa- 
tions coupled to a perfect fluid and a fc-essence field on a time interval 
[to,ti) with ti < oo. Suppose that the fc-essence field satisfies the dominant 
energy condition. If H{t) is bounded and if (0,0^/2) remains in a compact 
subset of G then the solution can be extended to a longer time interval. 
Proof The proof follows the discussion in section 4 of [20] . The equations of 
motion of the fluid and the fc-essence field can be written in the form 

A{z, gij, kij)dz/dt = F{z, gij, kij) (26) 

where z = (p, m\ 0, 0) and tensors are expressed in components with respect 
to a left invariant frame of the Lie algebra defining a given Bianchi type. The 
matrix A is invertible. The standard local existence theorem for ordinary 
differential equations [12] applies to this system. It gives local existence 
of solutions and shows that the solution can be extended to a longer time 
interval provided {z,gij, kij) remains within a compact subset of the domain 
of definition of the coefficients. It will be shown that this follows from the 
assumptions of the Lemma. A first step is to show that under the given 
assumptions gij, {detg)~^ and kij remain bounded. The proof adapts an 
argument used in the proof of Lemma 2.1 of jTH]. There the boundedness 
of the integral of kijk"^^ was deduced from certain energy conditions. In fact 
the dominant energy condition alone is enough. Equation (jSI) can be used to 
bound the integral of cxijcr*-' where aij is the tracefree part of kij. Together 
with the assumed boundedness of H this gives the desired conclusion. It 
remains to show that p, p~^ and are bounded.. This can be done just as 
in [2n|. 

As in |2ni! this result can immediately be extended to the case of several 
non-interacting fluids, e.g. a mixture of dust and radiation. It can also be 
adapted to the case of collisionless matter. For this purpose the existence 
theorem of section 2 of ^H] has to be extended but this is straightforward. 
The quantities and <p have to be added to those for which an iteration is 
carried out. The result is: 
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Lemma 2 Consider a solution of Bianchi type I- VIII of the Einstein equa- 
tions coupled to a collisionless gas described by the Vlasov equation and a 
/c-essence field on a time interval [to, ti) with ti < oo. Suppose that the 
fc-essence field satisfies the dominant energy condition. If H{t) is bounded 
and if (0, 0^/2) remains in a compact subset of G then the solution can be 
extended to a longer time interval. 
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